Long range scattering resonances in strong-field seeking states of polar molecules 
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We present first steps toward understanding the ultracold scattering properties of polar molecules 
in strong electric field-seeking states. We have found that the elastic cross section displays a quasi- 
regular set of potential resonances as a function of the electric field, which potentially offers intimate 
details about the inter-molecular interaction. We illustrate these resonances in a "toy" model 
composed of pure dipoles, and in more physically realistic systems. To analyze these resonances, 
we use a simple WKB approximation to the eigenphase, which proves both reasonably accurate and 
meaningful. A general treatment of the Stark effect and dipolar interactions is also presented. 
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I. INTRODUCTION 

Conventional spectroscopy of atoms and molecules be- 
gins with the premise that the energy levels of the species 
being probed are fixed although they may shift slightly 
in electromagnetic fields. These levels are then interro- 
gated by energy-dependent probes such as photons or 
charged or neutral particles. Information on the ener- 
getics and structure of the molecules is extracted from 
absorption energies, oscillator strengths, selection rules, 
etc. In these investigations, the study of resonances has 
played a central role. 

With the advent of ultracold environments for atoms 
and molecules, this general view of spectroscopy can be 
inverted. Cold collisions provide a nearly monochromatic 
probe of near-threshold intcrmolccular interactions, with 
resolution set by the milliKclvin or microKelvin tem- 
perature of the gas. In this case the energy levels of 
nearby resonant states can be tuned into resonance with 
the zero-energy collisions. For cold gases of alkali atoms, 
this strategy is already in widespread use. It exploits the 
fact that the Zeeman effect can shift the internal energies 
of the atoms over ranges orders of magnitude larger than 
the collision energy itself. In this way, atoms can be made 
to resonate when they would not naturally do so (i.e., in 
zero field) . Measurement of these "Feshbach" resonances 
(more properly, Fano-Feshbach resonances) has yielded 
the most accurate determination of alkali-alkali potential 
energy surfaces for near-threshold processes Q . 

The current experimental and theoretical push to cre- 
ate and study ultracold molecules Q will lead to many 
more opportunities for this novel kind of spectroscopy, 
since molecules possess many more internal degrees of 
freedom than do atoms. There will be, for example, nu- 
merous Fano-Feshbach resonances in which one or both of 
the collision partners becomes vibrationally or rotation- 
ally excited [lEIU; statistical arguments suggest that 
these resonances will be quite narrow in energy, a fact re- 
lated to their abundance |3|. A second class of resonant 
states will occur when the constituents are excited into 
higher- lying fine structure or hyperfine structure states, 
more reminiscent of the resonances observed in the alkali 
atoms. In molecules, these resonances are naturally also 



tunable in energy using magnetic fields jg- 

In this paper, we are primarily interested in a third 
class of resonance: potential resonances that are engen- 
dered by altering the intermolecular potential energy sur- 
face itself. This capability becomes especially prominent 
in cold collisions of heteronuclear polar molecules, whose 
dipolar interactions are quite strong on the scale of the 
low translational temperatures of the gas. At the same 
time, the dipole moments of the individual dipoles can 
be strengthened or weakened as well as aligned by an 
applied electric field. 

This kind of resonance includes the shape resonances, 
discussed in the context of cold atoms polarized by strong 
electric fields A second set, dubbed "field-linked" 

resonances, has been studied in some detail in Refs 
[H> EH lu| • These states appear in potential energy sur- 
faces (PES's) that correlate to weak electric field-seeking 
states of the free molecules. They are weakly bound, 
long-ranged in character, and indeed do not appear to 
exist without an electric field present. Similar resonances 
are predicted to occur in metastable states of the alkaline 
earth elements at low temperatures 113 • 

A rich set of potential resonances emerges among 
strong-field seeking states, and this is the subject of 
the present paper. Strong-field seekers are of increas- 
ing importance experimentally, since they would enable 
molecules to be trapped in their absolute ground states 
where no 2-body inelastic collision processes are available 
to harm the gas. In this case, the colliding molecules are 
free to approach to within a small internuclear distance of 
one another; the resulting potential resonances therefore 
can probe detailed intermolecular dynamics near thresh- 
old. The resulting data, consisting of scattering peaks as 
a function of electric field, can be thought of as a kind of 
"Stark spectroscopy" . Just such a tool has been applied 
previously in precision measurements of alkali Rydberg 
spectra 

In this paper we explore such spectra in ultracold 
molecules, finding that the spectrum is dominated by 
a quasi-regular series in the electric field values. Such 
a series is the fundamental building block of molecular 
Stark spectroscopy and plays a role analogous to the 
Rydberg series in atomic spectroscopy. In both cases, 
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the series lays out the fundamental structure of the un- 
perturbed, long-range physics between interacting enti- 
ties. In the case of atoms, the deviation from an un- 
perturbed Rydberg series, encapsulated in the quantum 
defects, yields information on the electron-core interac- 
tion 0, 0|. Similarly, it is expected that differences 
in observed Stark spectra from those presented here will 
probe the short-range intermolecular interactions. 

To emphasize the generality of these resonances, we 
first briefly develop in Sec. II a general formalism 
for arbitrary polar collision partners, covering explicitly 
Hund's cases (a) and (b) as well as asymmetric rotor 
molecules. Using this formalism, we introduce in Sec. 
Ill the potential resonances using a simplified version of 
the molecular gas in which all dipolcs arc assumed to 
be perfectly aligned and where molecular fine structure 
plays no role. In Sec. IV, we consider the case of more 
realistic molecules, where fine structure does intervene. 
We show that the structure of the potential resonances 
is unchanged, but that additional narrow Fano-Feshbach 
resonances do appear. Throughout the article we em- 
phasize how the various types of resonance can be classi- 
fied and organized by simple considerations involving the 
WKB approximation. 

II. GENERAL FORM OF THE STARK EFFECT 
AND DIPOLAR INTERACTION 

In this section, we briefly recapitulate our scattering 
model which is presented in Refs.[l(l[lftllfl]. The molec- 
ular scattering Hamiltonian can be written as 

H = T Q + V S R + V L r + H s + H int (1) 

Where To is the kinetic energy operator, Vsr is the short 
range potential energy surface (PES), Vlr is the long 
range interaction of the polar molecules, H$ is the Stark 
Hamiltonian, and Hint is the Hamiltonian describing the 
internal degrees of freedom of the molecules. For the 
purpose of this paper we disregard Vsr- The long range 
interaction is dominated by the dipole-dipole interac- 
tion [lOj. Furthermore an electric field can significantly 
change the structure of the interaction and be used to 
control the molecular collisions 0, 0] . 

The first step in performing a scattering calculation is 
to construct a basis of molecular energy eigenstates in an 
external field. Here we generalize the approach slightly 
to make it applicable to various types of polar molecules. 
Once the Stark Hamiltonian is obtained, the result can 
be used in the Hamiltonian of the dipolar interaction, 
showing their common physical origin in terms of electric 
fields. 

To simplify the analysis and calculations, we assume 
that the vibrational degrees of freedom are frozen out at 
low temperatures; hence can treat the molecules as rigid 
rotors. The molecular state is described in terms of the 
state | JMjfl), where J is the molecule's rotational plus 
electronic angular momenta, Mj is the projection of J 



onto the lab axis, and f2 is J's projection onto the molecu- 
lar axis. To describe the rigid rotor molecular wave func- 
tion, we use (a,P,j\JMjQ) = yJ^rD^ jSl {a, /3, 7), 
where a, /3, and 7 are the Euler angles defining the molec- 
ular axis and Df Ijn is a Wigner D-function |20| . 



A. The Stark effect 

Since the electric field is a true vector (as opposed to 
a pseudovector), it only couples states of opposite parity. 
This implies that the Stark energies vary quadratically 
with low electric field; they vary linearly only at higher 
fields once the Stark energy is greater than the energy 
splitting of the two states. The Stark Hamiltonian has 
the form 

H s = -$■£. (2) 

For the current discussion, we pursue a general approach 
and allow the field to point in any direction in the lab- 
oratory reference frame. This general approach allows 
the matrix elements to be used in the dipolar interac- 
tion. However when explicitly considering the molecular 
states in an electric field, we that the field lies in the z 
direction. 

The Stark interaction can be evaluated by decompos- 
ing the electric field in spherical coordinates and then 
rotating the dipole into the lab frame. In this repre- 
sentation, the Stark Hamiltonian has the form Hs = 
— ^ g /i£D^o, where fj, is the electric dipole moment of 
the molecule and D 1 * is a Wigner D-function. D]^ is 
equal to (— l) 2q C^, a reduced spherical harmonic |21| . 
and q is the projection quantum number of £ onto the lab 
axis. At the heart of evaluating the Stark effect, we find 
the operator D 1 * coupling two molecular states, which 
themselves are described by D functions. The Stark ma- 
trix element is an integral of three D-functions over the 
molecular coordinates, averaging over molecular orienta- 
tion. The evaluation of the integral results in selection 
rules whose details depend on the molecular specifics. 
For a valuable qualitative discussion on the Stark effect 
see [22^. 

A given molecule may also have a nuclear spin that 
generates a hyperfine structure. In this case, it is more 
appropriate to present the matrix elements in the hyper- 
fine basis, where F and Mp define the state. Here F is 
the sum of J and the nuclear spin / and Mp is the pro- 
jection of F in the lab frame. We use the Wigner-Eckart 
theorem to compute the Stark matrix elements in a com- 
pact form. The matrix elements of the Stark effect are 
written as 

(aFM F \H s \a'F'M F ) = -fi£(aFM F \Dl^\F' M' F a'), (3) 
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which contains a purely geometrical matrix element 



aFM F \D%\F'M' F a') = [F^-l) 1 



+M F +F' 



F' 

M' F M F 



M' F 



F 

-M f 



(aF\\Di*\\a>F'). (4) 



Here (aF \\ Dg* \\ a'Fl) is the reduced matrix element 
and a represents all remaining quantum numbers needed 
to uniquely determine the quantum state, [j] is a short- 
hand notation for y/2j + 1. We have left q in the matrix 
element simply as a place holder. Its value is assumed 
to be q = Mp — M' F in accordance with conservation of 
angular momentum. If we consider the case where the 
electric field points in the z direction, then q — 0, and we 
find that Mp = M' F , This means that Mp is a conserved 
quantum number. We have used the convention of Brink 
and Satchler to define the Wigner-Eckart theorem and 
reduced matrix elements l2lL 



where the index a represents the set of quantum numbers 
N and S. 

This formalism is easily extended to include asym- 
metric rotors by including the rotational Hamiltonian 
to construct molecular eigenstates. For an asymmetric 
rotor, there are three distinct moments of inertia and 
therefore three distinct rotational constants. The rota- 
tional Hamiltonian is H rot = A3 2 a + B3\ + CJ 2 , where 
A > B > C and a, b and c are the axis labels in the molec- 
ular frame. This additional structure mixes Q such that 
it is not a good quantum number, implying we that need 
to diagonalize the rotational Hamiltonian along with the 
Stark Hamiltonian to obtain the molecular eigenstates. 
We consider the primary effect of this additional struc- 
ture to change the progression of rotationally excited 
states. For a complete discussion on asymmetric rotors 
seeRefs. HHH. 



B. Molecular examples 

We present a few examples of reduced matrix ele- 
ments for specific molecular symmetries. First, consider 
a Hund's case (a) molecule with ft ^ 0. The OH radi- 
cal, with ground state 2 n 3 / 2 , is a good example of this. 
The energy eigenstates of such a molecule in zero elec- 
tric field are eigenstates of parity, i.e., \JMjCle) where 
e = +(—) represents the e (/) parity state and Ct = |fi|. 
The parity of this molecule is e(— 1) J_1 / 2 if J is a half 
inte ger or e(— 1) J if J is an integer. (For details see Refs. 
fl8t Il9t Taking this molecular structure into ac- 

count, we find that the reduced matrix element is 



(aF 



n 1 * 



F'a 1 ) = (-) 1+/ + F + J+J '- f2 [F', J, J'] 
F F' 1 1 / J'_ 1 J \ 

j' j i j \ -n o n ) 

'l + ee '(_l)W+2U+l\ 



(5) 



Here the index a represents the set of quantum num- 
bers e, fl and J. We have introduced another notation: 
\ji,h,~.,3N] = v / (2ji + l)(2j 2 + l)---(2j N + l). This 
reduced matrix element is the same for any case (a) 
molecule. 

On the other hand consider a Hund's case (b) molecule 
with L — 0. Many molecules fit this mold such as 
hctcronuclear alkali dimers and SrO with 25+1 E ground 
states. Here the parity of a state is directly identified by 
the value of J, where parity = (— 1) . The Stark effect 
therefore mixes the ground state with the next rotational 
state. We find the reduced matrix element to be 



(aF || D U || F'a') - ^f+i+J+J'+s+N+N' 

F F' 1 
J' J I 



x[N, N',J,J',F'} 



J J' 1 
N' N S 



N' 1 N 




(6) 



C. Dipole-dipole interaction 

The intrigue of polar molecules is their long range 
anisotropic scattering properties whose origin is the 
dipole-dipole interaction. The interaction in vector form 



3(R ■ Ai)(R ■ Ag) - £i ■ £2 
R 3 



(7) 



where pn is the electric dipole of molecule i, R is the in- 
termolecular separation, and R is the unit vector defining 
the intermolecular axis. This interaction is conveniently 
rewritten in terms of tensorial operators in the laboratory 
frame as: 



-—Ti- 



•1)«C* • (Mi® I*) 2 ,- 



(8) 



Here C^_ q (9, <j>) is a reduced spherical harmonic that acts 
on the relative angular coordinate of the molecules, while 
(fii ® is the second rank tensor formed from the two 
rank one operators determining the individual dipoles. 
For this reason, matrix elements of the interaction are 
also given conveniently in terms of the matrix elements in 
Eq. (0J. The matrix elements of the dipolar interaction 
are: 



(nirmlH^ll'^l'm't) 



f_^\— mi+m'j+l I L 1 



R :i 



x(lmi\ C( mi _ m , ) \l'm' l )(a 1 F 1 M Fl \ D'* \ a\F[M' Fi } 
x(a 2 F 2 Mp 2 \D^ \a' 2 F^M'p 2 ) 

1 M Fl - M' Fi Mp 2 - M' F2 mi -ml) (9) 



where 



(lm l \q mi _ m , ) \l'm' l ) = (-ir'[l,l'}[ l Q I [ 



I 



I' 



-mi mi — m\ m\ 



(10) 



Equation shows that once the Stark Hamiltonian 
has been constructed for a particular molecule, then the 
Hamiltonian describing the dipolar interaction can be 
achieved with little extra effort. This result reflects the 
fact that in Eq. each dipole is acted on by the elec- 
tric field of the other. We have also used a shorthand 
to represent the channel, |12Zto/), where |f) denotes the 
quantum state of the first molecule and as |2) for second 
molecule. At this point, we disregard all interactions be- 
tween molecules except the dipolar interaction. 

To exploit an analogy with Rydberg atoms, the long 
range Coulombic interaction is well understood. With 
this understanding of the long range physics, a solution 
to the complete problem is achieved by matching it to 
the short range solution or a parameterization of the 
short range interaction [l7j |. This idea was, in fact, im- 
plemented as a numerical tool in Ref.@ , which dealt with 
collisions of atoms in strong electric fields. To this end, 
we pursue an understanding of the long range character- 
istics of the dipolar scattering. We then envision merging 
the long range physics with the short range physics, or 
a parametrization thereof, to offer insight into the short 
range interaction and to explore the dynamic interaction 
of polar molecules. With this idea in mind we first ex- 
plore "pure" dipolar scattering. 



III. DIPOLAR SCATTERING 

Our primary interest in polar molecule scattering is 
how the strong anisotropic interaction affects the sys- 
tem. As a first step illustrating the influence of dipolar 
interactions, we present a simple model composed of po- 
larized dipoles with no internal structure. Strictly speak- 
ing, this system is created by an infinitely strong electric 
field that completely polarizes the molecules and raises 
all other internal states to experimentally unattainable 
energies. Thus the only label required for a channel is 
its partial wave, I = {0,2,4,...} in the numerical ex- 
ample given here. The matrix elements of the dipole- 
dipole interaction are taken to be (12/0|if MM |12Z'0) = 
— 0.32/i 2 (Z0|CQ \1'0), which is typical for molecules like 
RbCs or SrO. We then artificially vary the dipole mo- 
ment fj, to see the effect of an increasingly strong dipolar 
interaction. Pragmatically speaking, varying /i parallels 
changing the electric field. The intention of this model 
is to focus on the effect of direct anisotropic couplings 
between the degenerate channels, as measured by their 
effect on the partial wave channels. 

For this model, we use a reduced mass of m r = 10 4 
a.u., typical of very light molecules. We moreover as- 
sume that the molecules approach one another along the 
laboratory z-axis, so that only the projection to; = 
of orbital angular momentum is relevant. To set a con- 
crete boundary condition at small R, we apply "hard 
wall" boundary conditions, ip(Ri n ) = at a character- 
istic radius Ri n = 20 a.u. (In Sec. V we will relax this 
restriction.) We pick the collision energy to be nearly 
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FIG. 1: (a) Adiabatic curves of the pure dipole system for 
two values of the dipole moment /j,: y, = 0.1 (gray) and /i = 1 
(black) (b) The cross section of the polarized dipole model 
versus fi, at a collision energy E = 10" 12 K. The brackets 
denote predicted resonant positions using the adiabatic WKB 
phase (AWP) approximation. The black bracket represents 
for the phase contribution from the lowest adiabatic curve 
and the gray is from all remaining contributions. 



zero, namely 10 -12 K. To converge the calculation for this 
model requires inclusion of partial waves up to I = 14, 
and numerical integration of the Schrodinger equation 
out to R — Rao = 1 x 10 5 a.u., using the log-derivative 
propagator method of Johnson [25| . 

To get a sense of the influence of increasing the dipole 
moment, we first look at adiabatic curves of the system. 
Figure n (a) shows two different sets of adiabatic curves: 
a gray set with /i = 0.1 (a.u.) and a black set fi = 1.0 
(a.u.). In each set, the four lowest adiabatic curves are 
shown. Looking at these curves, we can see two charac- 
teristic effects of increasing fi. First, the lowest curve be- 
comes much deeper. Second, the higher adiabatic curves, 
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originating form non-zero partial waves, may support 
bound states at short distance states, i.e. within the cen- 
trifugal barrier. Both these effect may generate bound 
states, leading to distinct classes of scattering resonances 
as ji is varied. The deepening of the lowest adiabatic 
curves induces potential resonances, whereas the higher- 
lying curves lead to narrow shape resonances, wherein the 
molecules must tunnel through the centrifugal barrier. 

The different classes of resonances can clearly be seen 
in cross section, as shown in Fig. ^ (b). The broad 
quasi-regular set of resonances seen in the cross section 
are the potential resonances originating from the lowest 
adiabatic curve. The narrow shape resonances appearing 
sporadically in the spectrum originate from the higher- 
lying curves. For the purpose of this paper, we focus on 
the wide potential resonances and simply acknowledge 
the existence of the narrow shape resonances. 

To show that the broad resonances primarily belong 
to the lowest potential and the narrow shape resonances 
belong to the higher-lying potentials, we use an eigen- 
phase analysis. The eigenphase can be thought of as the 
sum of the phase shifts for all of the channels; thus it 
tracks the behavior of all the channels simultaneously. 
The eigenphase is defined as 

4> e i gen = ^tan-\\f). (ff) 

i 

Here Xf are the eigenvalues of the K matrix from the full- 
scattering calculation j2(|. (The K matrix is related to 
the more familiar scattering matrix by S = (1 + iK) /(l — 
iK).) When the system gains a bound state, it appears 
as a 7r jump in eigenphase. The eigenphase of the system 
is shown in Fig. 121(a), as the solid line with many abrupt 
steps. 

To analyze this situation further, we construct an ap- 
proximate eigenphase as follows. First, we asses the to- 
tal phase accumulated in each adiabatic channel, using a 
WKB prescription: 

^wkb(p) = y/-2m r V^(ti,R)/^. 

J Rin 

Here (i) stands for the i th adiabatic curve. For the low- 
est adiabatic curve, which is always attractive, the range 
of integration is from R in to R out = oo. For high- lying 
channels that possess a barrier to scattering at zero col- 
lision energy, the limits of integration are from R in to 
the inner classical turning point of the barrier. This will 
yield some information on shape resonances trapped be- 
hind the barrier, but we will not make much of this in the 
analysis to follow. Finally, we add together the individ- 
ual WKB phases to produce an approximate eigenphase 
shift, which we dub the "adiabatic WKB phase" (AWP): 

4>wkb{^) = ^Wk-bW- ( 12 ) 




FIG. 2: (a) Eigenphase (solid) and total adiabatic WKB 
phase (AWP,dashed) for the dipole-scattering model, (b) 
AWP contributions for the lowest adiabatic curves (4>wkb) 
with the black bracket indicating when <Pwkb P asses through 
a multiple of -n. The two remaining curves are contributions 
with non-zero partial waves {4>w'kb)- Phe S ra y bracket in- 
dicates where the sum of these contributions pass through a 
multiple of 71". 

Since we are not concerned with properties of the phase 
associated with the higher-lying adiabatic curves, we do 
not consider the connection formula now. 

The total AWP for this system is shown as a dashed 
line in Fig. [21 (a) . It tracks the eigenphase well but of- 
fers more information if we decompose the AWP into its 
contributions. Figure [21 (b) shows the individual con- 
tributions of the sum. The largest contribution is the 
phase from the lowest adiabatic curve, which can be as- 
sociated with potential resonances. A black bracket ap- 
pears above this phase contribution with vertical marks 
indicating when it passes through an integer multiple of 
7r, i.e., when we expect to see a potential resonance in 
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the cross section. This same bracket is plotted in Fig. Q] 
and shows good agreement between the locations of the 
potential resonances and the AWP predictions. We con- 
clude from this that the main resonance features in the 
Stark spectrum arise primarily from this single potential 
curve. 

In Fig. [21 (b) the two remaining phase contributions 
originate from non-zero partial wave channels that pos- 
sess centrifugal barriers; see Fig. db) . The gray bracket 
indicates where the sum of these two contributions pass 
through an integer multiple of tt, and thus represent a 
guess for where the shape resonances lie. This gray 
bracket is also shown in Fig. ^ The agreement with 
the position of the narrow resonance features in the cross 
section is not nearly as good as it is with the broad poten- 
tial resonances. This indicates a more involved criterion 
for shape resonances. Nevertheless, the AWP predicts 
15 shape resonances, and there are 13 in the range of /i 
shown. The AWP appears to offer a means to roughly 
predict the number of shape resonances in this system, 
even though it does not predict the locations exactly. 

A main point of this analysis is that the AWP in the 
lowest adiabatic channel alone is sufficient to locate the 
potential (as opposed to shape) resonances, without fur- 
ther modification. For the rest of this paper, we focus 
on the potential resonances in more realistic molecules 
with internal molecular structure. The general analysis 
in terms of a single-channel AWP will still hold, but an 
additional phase shift will be required to describe the 
spectrum. 



IV. STRONG FIELD SEEKERS 

Strong-field-seeking molecules are approximately de- 
scribed as polarized in the sense of the last section, be- 
cause their dipole moments are aligned with the field. 
They will, however, contain a richer resonance structure 
owing to the presence of low-lying excited states that can 
alter the dipolar potential energy surface at small- R. 

For concreteness, we focus here on molecules with a 
1 E ground state. Heteronulcear alkalis fit into this cat- 
egory and are rapidly approaching ground state produc- 
tion with various species |27|. As examples, we pick RbCs 
and SrO in their ground states. Ground state RbCs has 
been produced experimentally |28j| . As for SrO, promis- 
ing new techniques should lead to experimental results 
soon [2^. For simplicity we include only the J = 0, J = 1 
rotational states, and freeze the projection of molecular 
angular momentum to Mj = 0. This restricts the num- 
ber of scattering thresholds to three, identified by the 
parity quantum number of the molecules in zero field. 
The parity quantum numbers for the three thresholds 
are ( , — +, ++}• This model is similar to the one pre- 
sented in Ref. [l(| which can be easily constructed for 
any rigid rotor when only including two molecular states. 
One immediate consequence of multiple thresholds is the 
presence of rotational Fano-Feshbach resonances in the 




R(a.u.) 



FIG. 3: (a) Adiabatic curves for the RbCs system at two dif- 
ferent field values, £ = 0(V/cm) (gray) and £ — 5000(V / cm) 
(black) . (b) Log-log plot of the absolute value of the two low- 
est adiabatic curves The blue dashed line is proportional to 
l/R 3 . 



collisional spectrum. 

The first example is RbCs, whose physical parameters 
are /z = 1.3 D, m r = 110 a.m.u., and B e = 0.0245(if) 
|28|. As before, we apply vanishing boundary condition 
at Ri n — 20 a.u. To converge this calculation over the 
field range considered, we require partial waves up to 
I = 30. We first look at the adiabatic curves of the system 
to get an understanding of how the real system deviates 
from the simple model presented above. In Fig. 0(a) we 
plot the 6 lowest adiabatic curves for the RbCs system 
with only four partial waves, so the figure is more easily 
interpreted. The sets of adiabatic curves shown are for 
two different field values: the gray set has £=0 and the 
black set has £ = 5000 (V/cm). There are two important 
features that differ from the dipole example. First, there 
are two higher thresholds, and the electric field shifts 
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these apart in energy as the field is increased. Second, the 
electric field dramatically changes the radial dependence 
of the Hamiltonian. 

The difference in thresholds can be seen clearly in Fig. 
13 (a) where the lowest excited threshold moves from 
0.05(10 at £=° t0 0.18(10 at £=5000 (V/cm). The 
difference in radial dependences for the two cases is seen 
more clearly in a log-log plot of the two lowest adiabatic 
curves for both fields, as shown in Fig. [3] (b). The gray 
set corresponding to zero field has two distinct asymp- 
totic radial powers. At large R the lowest adiabatic curve 
has a 1 /R e behavior asymptotically because of couplings 
with channels far away in energy [~ 0.05(10]. However 
as R approaches zero the dipolar interaction has over- 
whelmed the rotational energy separation and the radial 
dependence becomes 1/R 3 in character at about R=100 
(a.u.). For reference, the dashed line is proportional to 
1/R 3 . With £= 5000 (V/cm), the two black curves show 
the radial dependence of the adiabatic curves. The lowest 
curve now has nearly 1/R 3 over the whole range shown. 
Asymptotically when the centrifugal barrier is larger, the 
radial dependence will change to 1/R 4 ^(j- The second 
lowest adiabatic curve is also significantly altered by the 
strong dipolar interaction as can be seen in Fig. |3| (b) . 

In Fig. 0] (a) , we plot the cross section for the model 
RbCs system. This spectrum is riddled with narrow 
Fano-Feshbach resonances, but is still dominated by a 
series of potential resonances similar to the one in Fig. 
n (b) . There are two sets of AWP predictions shown as 
over-brackets. To understand their difference, we look 
to Fig. E|(b). The AWP for the lowest adiabatic curve 
is shown in Fig. 0] (b) for two different zero-field phase 
values. The gray curve is the AWP that is directly com- 
puted from the method described above, Eq. 1121) . The 
locations where it passes through an integer multiple of 
7T are indicated by the gray triangles. Referring back to 
Fig. Q] (a), where the same gray triangles appear we see 
that this simple estimate does not reproduce the reso- 
nance position. 

We can, however, introduce an additional overall phase 
shift to account for the difference in short-range interac- 
tions from the pure polarized case. The shifted AWP 
reads 



0wkb(£) 



TfSdefe 



(13) 



By treating Sdefect as a fitting parameter, we can obtain 
the resonance positions indicated by the black bracket 
in Fig. 01 which agree quite well with the resonance 
positions in the close-coupled calculation. To do so re- 
quires, in this case, a phase shift Sdefect = 0.14. In anal- 
ogy with Rydberg spectroscopy we consider the shift we 
have added to be a "quantum defect" that accounts for 
the effect of the short-range interaction. The additional 
phase shift reflects the influence of short range physics on 
the scattering, such as curve crossings with curves from 
higher thresholds. 

The AWP also saturates with field, as can be seen in 
Fig. 0| (b) . This occurs because the electric field eventu- 




' 1000 ' 2000 3000 ' 4000 ' 5000 
Electric Field (V/cm) 




1000 2000 3000 4000 5000 
Electric Field (V/cm) 



FIG. 4: (a) Cross section for RbCs in strong-field seek- 
ing states, including the 10 potential resonances that the 
AWP predicts (solid bracket). The AWP phases are shown 
in (b). The black bracket corresponds to the AWP pre- 
dicted resonances for the lowest adiabatic curve after a defect, 
Sdefect = 0.14, has been added according to Eq. 1131 . The 
gray triangles are the AWP predicted resonances without the 
defect. 



ally fully polarizes the molecules, so the dipole moment 
cannot increase further. The effect can also be seen in 
the spacing of the potential resonances. At low fields, 
the potential resonances occur frequently in field. Then, 
as the field is further increased, the resonances occur less 
often in field, which is a signature of dipole moment satu- 
rating and therefore an increasing field having less effect 
on the molecular interaction. 

As a second example, we consider SrO, which has the 
physical parameters m r =52 a.m.u, /Lt=8.9 D, and B e = 
0.5 K [23. We choose this molecule for its compara- 
tively large mass and dipole moment, which guarantee a 
large number of resonances. Fig. [SJwe have plotted the 
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1000 2000 3000 4000 5000 
Electric Field (V/cm) 



FIG. 5: Cross section for SrO with a bracket indicating res- 
onance positions predicted by the phase shifted AWP for the 
lowest adiabatic curve. There are 33 PR shown in this range 
of electric field. Sdefect = 0.215 

cross section for SrO, which is dominated by the quasi- 
regular potential resonance series. As before, the black 
bracket indicates where the phase shifted AWP predicts 
the potential resonances, and we see the agreement is 
good. Furthermore, the series has not terminated since 
we have not completely polarized the dipole. The series 
of potential resonances saturates at 17.5 (kV/cm) after 
a total of 43 potential resonances have been induced. To 
line up the AWP's predictions and the actual cross sec- 
tion requires a defect of Sdefect — 0.215. 

We have picked two examples to illustrate how the po- 
tential resonances will appear in the context of collisional 
spectroscopy. These resonance will occur to varying de- 
grees in the strong-field seeking collisions of all polar 
molecules. For example, we can also make similar predic- 
tions for an asymmetric rotor molecule such as formalde- 
hyde (H2CO). We find that this molecule should possess 
six potential resonances in the field interval from to 50 
(kV/cm). 

It is worth noting that portions of similar resonance 
series was anticipated in cold atomic gases subjected to 
electric fields 0- However very few such resonances are 
likely to be observed, owning to the enormous fields (~ 
MV/cm) required to generate them. In polar molecules, 
by contrast, the entire series should be readily observable. 

V. COLLISIONAL SPECTROSCOPY 

Through the course of this work we have shown 
that the zero-energy cross section of weak-field-seeking 
molecules is dominated by a set of broad potential reso- 
nances. Even though these resonances are themselves in- 
triguing, their properties can be exploited to learn much 
more about the system. The general structure of the po- 
tential resonances is governed by the long range dipolar 
interaction, which has a predictable and common form. 



With a clear understanding of this interaction and how it 
induces resonances, it could be exploited to learn about 
the short range interaction of the molecules. This is be- 
cause details of where the lines appear must also depend 
on the boundary condition experienced by the wave func- 
tion at small values of R. Therefore the spectrum con- 
tains information on the small R intermolecular dynam- 
ics. Thus by studying the potential resonance series we 
can extract information about the short range dynamics. 

This idea is similar to quantum defect theory, which, 
for example, has been very successful in Rydberg spec- 
troscopy. The short range physics of the electron inter- 
acting with the nucleus is complicated and not easily 
solved. However once the electron is out of the small 
R region, it enters into a pure coulomb potential where 
the motion of the electron is well understood. The effect 
of the short range must be merged with the long range 
physics to form a complete solution. To account for the 
short range interaction, the energy can be parameterized 
by replacing the principal quantum number by an effec- 
tive quantum number n* = n — fx. This procedure is 
tantamount to identifying an additional phase shift due 
to the interaction of the electron with the atomic core. 
The idea of merging standard long range physics with 
complicated short range behavior has been applied suc- 
cessfully not only in Rydberg states of atom s Il7l l30| and 
molecules l31l . but also in atomic collisions |32j , cold col- 
lisions |33l l34l l35| , and dipole-dipole interactions of the 
type we envision here 

As a simple expression of this idea, we can alter the 
boundary condition applied at i?; n when performing the 
full scattering calculation and note its influence on the 
field dependent spectrum. For the above calculations 
we have imposed the standard vanishing boundary con- 
dition, ip(Ri n ) — for all channels. We now replace 
this condition with a uniform logarithmic derivative, 
b = (j^iJj)^^ 1 , at Ri n . Thus previously we set b = 00, 
but now we allow b to vary. The log-derivative is conve- 
niently represented as a phase: 

b = cotan(7r/3), (14) 

where p can lie between zero and one, covering all values 
of b from —00 to +00. For /3 = 0, the boundary condition 
is the one employed above, tjj{R in ) = 0, whereas for /3 = 
0.5, the boundary condition is ^V(^in) = 0. 

We have re-computed the collisional spectrum of RbCs 
for several different initial conditions, and plotted the 
field values of the potential resonances, E {b \ in Fig. © as 
sets of points. The values of (3 for the different calcula- 
tions arc (filled circle), 0.11 (filled square), 0.22 (filled 
triangle), 0.56 (hollow circle), 0.78 (hollow square), and 
0.89 (hollow triangle). The filled circles are resonant lo- 
cations for cross section in Fig. 0] (a). 

We next wish to demonstrate that each such spectrum 
can be identified by a single quantum defect parameter, 
as was done in the previous section. This entails picking 
a value of Sdefect and then setting the phase 4>wkb{^) 
equal to an integer multiple of 7r. This yields a set of 
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FIG. 6: The vertical curves represent electric field values 
where the AWP predicts potential resonances for a given de- 
fect, £wKB(Sdefect)- The points represent resonant field val- 
ues, in the full calculations with different initial bound- 
ary conditions. The values of b are given by cotan(7r/3) and (5 
for the different full numerical calculations are (filled circle), 
0.1 1 (filled square), 0.22 (filled triangle), 0.56 (hollow circle), 
0.78 (hollow square), and 0.89 (hollow triangle). Each set of 
values £ ™ is plotted at a height corresponding to its best-fit 

value Of Sdefect- 



With Sdefect = 0.14. 

We can compare the the resonant field values predicted 
by the AWP, Swkb (Sdefect), and resonant field values 
given by the full calculation with different initial condi- 
tions, E^ b \ To plot we have varied the height at 
which the set of points £^ is plotted until it aligns with 
£wkb{& defect)- Doing this we are able to to see how 
b and Sdefect are related. Thus in fig. we can see 
that even with different boundary conditions, the single 
AWP curve in fig 0] (b) can be used to predict the spac- 
ing between the potential resonances by varying a single 
parameter, Sdefect- This shows the AWP represents the 
long range scattering physics well, and that empirically 
extracted parameters like Sdefect will carry information 
about the short-range physics such as that embodied in 
b. 

VI. CONCLUSION 

A number of resonant processes may occur when two 
polar molecules meet in an ultracold gas. We have 
focused here on the dominant, quasi-regular, series of 
potential resonances between weak-field seeking states. 
These potential resonances originate in the direct defor- 
mation of the potential energy surface upon which the 
molecules scatter. Observation of these resonances may 
offer a direct means for probing the short range interac- 
tion between molecules. We have provided a means of 
analyzing this system with an adiabatic WKB phase in- 
tegral. This method shows how the system evolves with 
the application of an electric field. 



resonant field values, Zwkb{S defect)- Each Sdefect cor- 
responds to a particular approximate spectrum. The 
set of curves, £wkb(S defect), generated by continuously 
varying Sdefect are shown are shown in Fig. [H] as solid 
lines. The bracket in Fig. 0] (a) corresponds to the set 
of points where a vertical line intersects £wkb{S 'defect) 
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